NATIONAL  BUREAU  OF  STANDARDS 
microcopy  resolution  test  chart 


THE  DYNAMIC  LOOP  EFFECT  ON  THE 
MISSISSIPPI  RIVER.  PROJECT  DESIGN 
FLOOD  FLOW  LINE 


D D C 


Tho  finding*  in  dtl*  report  nr*  not  to  bn  constrvod  os  on  officiol 
DOpurtmsnt  of  A#  Army  position  unlost  so  dosignotod 
by  efwr  owdwrizod  dotwMwti. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


THE  DYNAMIC  LOOP  EFFECT  ON  THE  MISSISSIPPI  RIVER 
'PROJECT  DESIGN  FLOOD  FLOW  LINE  — - Z J 


Final  re' 


Carl  J/Huval 


».  PERFORMING  ORGANIZATION  NAME  AND  ADORESS 

U.  S.  Army  Engineer  Waterways  Experiment  Station 

Hydraulics  Laboratory 

P.  0.  Box  631,  Vicksburg,  Miss.  39180 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


U.  S.  Army  Engineer  Division,  Lower  Mississippi 
Valley 

P.  0.  Box  80,  Vicksburg,  Miss.  39180 


18.  SECURITY  CLASS,  (of  thla  rmport) 


Hlarant  from  Controlling  Off  lea) 


Unclassified 


ISa.  OECLASSI FICATION/ DOWN  GRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


IB.  SUPPLEMENTARY  NOTES 


19.  KCY  WORDS  (Contlnua  on  rararaa  a I da  If  nacaaaary  and  Idantlty  by  block  numbar) 


Mathematical  models 
Peak  floods 
Mississippi  River 
Unsteady  flow 


Design  flood 
Dynamic  loop  effect 
Flood  control 
Flow  characteristics 


M.  ABSTRACT  fC Mka  a mwtwm  aMa  IF  imlRR'  md  Mmmtltr  k y Mock  m—karj 

— ^ The  unusually  high  1973  and  1975  flood  seasons  on  the  Mississippi  River 
revealed  apparent  significant  reductions  in  flood-carrying  capacity.  In  addi- 
tion, a heretofore  unappreciated  phenomenon  known  as  loop  effect  appeared  to 
have  a large  impact  on  stage-discharge  relations.  As  a part  of  a general  re- 
analysis of  the  design  flood  flow  line,  a study  was  initiated  to  help  quantify 
the  dynamic  loop  effect.  The  purpose  of  the  study  was  to  estimate  the  loop 
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magnitude  at  the  design  flood  to  be  added  to  the  separately  evaluated  steady- 
state  design  flood  flow  line.  ^ 

A previously  developed  numerical  model  for  calculations  of  variable  stage- 
discharge  relations  was  used.  The  study  showed  that  this  model  could  not 
account  for  all  of  the  loop  magnitude,  probably  due  to  flow  complexity  and  mode] 
assumptions.  An  ad  hoc  method  to  simulate  the  total  loop  was  implemented  in- 
volving double-valued,  variable  Manning's  n corresponding  to  the  rising  or 
falling  phase  of  the  flow  hydrographs.  The  method  of  simulation  suggests  a 
higher  channel  resistance  at  falling  river  stages.  Extrapolation  to  the  design 
flood  is  based  on  the  assumption  of  similarity  of  river  behavior  to  the  1973 
flood.  This  assumption  is  consistent  with  the  steady  flow  line  computations 
using  river  resistance  extrapolation  from  the  1973  flood  data  to  the  design 
flood. 

The  modified  model  was  applied  to  three  representative  gaging  stations  to 
quantify  the  loop  magnitude  along  the  lower  river.  Results  of  the  study  indi- 
cate that  the  dynamic  loop  allowance  should  be  as  follows: 

Vicksburg  - 0.9  ft 
Helena  - 1.0  ft 

Baton  Rouge  - 0.5  ft 

In  summary,  an  allowance  of  approximately  1.0  ft  should  be  added  to  the  flow 
line  from  about  Helena  to  Red  River  Landing.  An  amount  of  0.5  ft  should  be 
sufficient  from  Red  River  Landing  to  about  Baton  Rouge. 
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PREFACE 


This  investigation  was  sponsored  by  the  U.  S.  Army  Engineer  Divi- 
sion, Lower  Mississippi  Valley  (LMVD),  and  funded  jointly  by  the  U.  S. 
Army  Engineer  Districts,  Memphis,  Vicksburg,  and  New  Orleans.  The  study 
was  conducted  during  the  period  March  1976  to  October  1978  by  the 
Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  under  the  general  supervision  of  Messrs.  H.  B.  Simmons, 
Chief  of  the  Hydraulics  Laboratory,  and  M.  B.  Boyd,  Chief  of  the 
Hydraulic  Analysis  Division.  A preliminary  report  presenting  final 
computational  results  at  the  Vicksburg  and  Helena  gaging  stations  and 
tentative  results  at  Baton  Rouge  was  presented  to  LMVD  in  May  1977. 

The  study  was  conducted  and  this  report  was  prepared  by  Mr.  Carl 
Huval.  Mr.  M.  B.  Boyd  contributed  substantial  technical  input  to  the 
study.  Special  assistance  in  providing  data  was  contributed  by 
Mr.  Malcolm  Dove,  LMVD,  as  well  as  personnel  from  the  Vicksburg,  Memphis, 
and  New  Orleans  Districts. 

Commander  and  Director  of  WES  during  the  investigation  and  the 
preparation  and  publication  of  this  report  was  COL  John  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 


To  Obtain 


cubic  feet  per  second 

0.02831685 

cubic  metres  per 
second 

feet 

0.3048 

metres 

feet  per  second 

0.3048 

metres  per  second 

feet  per  second  per 
second 

0.3048 

metres  per  second 
per  second 

miles  (U.  S.  statute) 

1.609344 

kilometres 

square  feet 

0.09290304 

square  metres 

THE  DYNAMIC  LOOP  EFFECT  ON  THE  MISSISSIPPI  RIVER 
PROJECT  DESIGN  FLOOD  FLOW  LINE 


PART  I : INTRODUCTION 


1.  Recent  large-scale  flooding  during  the  1973  and  1975  flood  sea- 
sons on  the  Mississippi  River  and  its  tributaries  indicated  that  appar- 
ently significant  reductions  had  occurred  in  the  flood-carrying  capacity 
of  the  river  in  some  reaches.1  These  changes  resulted  in  higher  stages 
for  given  discharges  and  could  have  a major  impact  on  the  degree  of 
flood  protection  provided  by  engineering  works  such  as  levees.  As  a 
result  of  the  1973  flood,  a reanalysis  of  the  project  design  flood  flow 
line  was  initiated  by  the  U.  S.  Army  Engineer  Division,  Lower  Missis- 
sippi Valley  (LMVD),  and  the  U.  S.  Army  Engineer  Districts,  New  Orleans, 
Vicksburg,  and  Memphis.  The  main  effort  of  the  district  studies  in- 
volved the  use  of  a steady-state  flow  model,  HEC-2,  Water  Surface  Pro- 

2 

files,  to  recompute  the  project  flood  flow  line  along  the  river.  The 
1973  flood  peak  stage  and  discharge  data  and  recent  channel  cross- 
section  measurements  were  used  to  recalibrate  and  then  extrapolate  to 
design  flows. 

2.  In  addition  to  apparent  reduction  in  channel  capacity,  dis- 
charge measurements  indicated  that  some  increase  in  stage  was  attribut- 

3 

able  to  loop  effect.  It  is  well  known  that  changing  discharge  during 
a flood  can  produce  a so-called  "hysteresis  loop"  in  the  stage-discharge 
rating  curve  such  as  the  one  shown  in  Figure  1.  For  this  type  of  rating 
curve,  there  can  be  several  different  stages  for  a given  discharge.  The 
lower  stage  value  is  associated  with  the  rising  limb  of  the  discharge 
hydrograph  and  the  higher  value  occurs  during  the  recession  of  the  dis- 
charge. A hysteresis  loop,  which  is  caused  by  a variable  energy  slope 
due  to  changing  discharge,  has  also  been  referred  to  as  a "dynamic  loop" 
since  it  occurs  because  of  the  changing  or  dynamic  nature  of  the  flood 
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DISCHARGE 


discharge.  A plot  of  several  recent  loop  rating  curves  at  Vicksburg  is 
shown  in  Figure  2.  Prior  to  the  1973  flood,  the  1950  flood  was  the  most 
recent  major  flood  event  on  the  Mississippi  River.  As  shown  in  Fig- 
ure 2,  an  adopted  rating  curve  had  been  used  for  design  based  on  the 
1950  flood.  The  1973,  197^,  and  1975  flood  rating  curves  are  signifi- 
cantly higher  than  the  one  for  1950.  As  shown  in  Figure  2,  the  stage- 
discharge  relation  can  be  very  complex,  since  it  is  dependent  on  several 
river  and  inflow  factors. 

3.  The  1973  flood  showed  that  the  project  design  flood  flow  line 
had  to  account  for  the  looping  effect  on  the  Mississippi  River.  For 
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Figure  2.  Vicksburg  recent  floods,  measured  rating  curves 


preliminary  design  purposes,  a dynamic  loop  allowance  of  about  1.7  ft* 
had  been  added  to  the  1973  design  flow  line  at  Vicksburg,  Miss. , and  an 
amount  of  2.1+  ft  at  Arkansas  City,  Ark.  Smaller  allowances  were  made 
at  other  gaging  stations  along  the  river.  The  provisional  allowances 
were  made  on  the  basis  of  field-measured  loop  magnitudes  during  the 
1973  flood  and  a study  of  past  rating  curves  along  the  river.  A flow 
line  review  study  requirement  by  the  LMVD  was  to  use  an  analytical  and 
hydrologically  consistent  method  for  computation  of  the  dynamic  loop 
effect. 

I 

Purpose  and  Scope  of  Study 

1*.  The  purpose  of  this  study  is  to  identify  dynamic  loop  and  other 
unsteady  flow  effects  during  the  1973  flood  and  to  estimate  the  loop 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  presented  on  page  3. 
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magnitude  at  the  design  flood.  The  loop  magnitude  is  the  amount  to  be 
added  to  the  steady-state  project  design  flood  flow  line  computed  using 
HEC-2  at  the  design  flood.  As  a part  of  the  study,  the  model  developed 
by  "read  was  implemented  for  use  in  providing  a consistent  and  improved 
calculation  method  for  determining  dynamic  loop  effect  for  adjusting  the 
steady-state  flow  line.  The  1973  flood  year  was  used  as  the  basic  cali- 
bration and  Vicksburg,  Miss.,  Helena,  Ark.,  and  Baton  Rouge,  La.,  gaging 
stations  were  investigated.  The  model  was  calibrated  against  the  1973 
flood  at  each  of  those  stations  and  computations  were  made  for  the  proj- 
ect design  flood  58A-EN.  In  order  to  determine  the  effects  of  consecu- 
tive loops  occurring  on  the  rising  limb  of  the  flood  hydrograph,  several 
discharge  hydrographs  were  run  to  simulate  other  potential  flood  hydro- 
graphs. The  two  hydrographs  described  below  were  selected  as  the  main 
input  flows  for  the  loop  computations:  (a)  combined  hydrograph  using 
the  top  part  of  the  hypothetical  project  flood  hydrograph  58A-EN  with 
the  observed  rising  part  of  the  1973  flood  hydrograph  and  (b)  the  hypo- 
thetical project  flood  hydrograph  58A-EN  alone,  which  has  one  single 
peak. 
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PART  II:  THE  MODEL 


Background 

5.  The  river  engineer  is  frequently  concerned  with  the  conversion 
of  flood  stages  at  a given  location  along  a stream  channel  into  corre- 
sponding discharges  or  vice  versa.  This  is  accomplished  by  a relation 
between  stage  and  discharge  which  applies  to  that  particular  gaging 
station.  Such  a relation  or  "rating  curve"  is  usually  developed  empiri- 
cally from  a number  of  previous  streamflow  measurements  and  the  corre- 
sponding stages.  Unfortunately,  the  observed  measurements  of  stage  and 
discharge  will  not  usually  form  a unique  relation,  i.e.,  a single  value 
of  stage  does  not  correspond  to  a single  value  of  discharge.  Deviations 
of  the  measurements  from  a single-valued  rating  curve  can  be  the  result 
of  such  factors  as: 

a.  Discharge  measurement  errors. 

b.  Shifting  control  due  to  scour  or  fill  and  alluvial  bed  form 
changes. 

c_.  Variable  water-surface  slope  along  the  channel  due  to  un- 
steady flows. 

d.  Local  effects  causing  variable  water-surface  slope  due  to 
backwater  from  a tributary  or  overbank  storage. 

£.  Water  temperature  effect  on  alluvial  bed  movement. 

6.  While  the  effects  due  to  errors,  sediment  movement,  and  tribu- 
tary and  overbank  flows  are  not  easily  quantified,  flow  accelerations 

due  to  unsteady  nonuniform  flow  can  be  computed.  In  a recent  report, 

1* 

Fread  presented  a relatively  simple  mathematical  model  that  retains 
most  of  the  dynamic  loop  effects.  The  method  was  proved  accurate  for 
calculating  the  dynamic  stage-discharge  relation  at  a number  of  river 
gaging  stations. 


Model  Details 


7.  A thorough  description  of  the  model 

It 

documentation  has  been  presented  by  Fread. 


with  the  computer  program 
The  basic  concept  is  to 
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model  the  stage-discharge  relation  for  a particular  gaging  station  loca- 
tion on  the  river,  i.e.,  equation  terms  involving  spatial  gradients  of 
discharge  and  flow  depth  are  transformed  into  time  variations  at  the 
location.  The  model  development  requires  several  assumptions  that  are 
applicable  for  a short  section  of  channel  containing  the  gaging  station 
or  forecast  point: 

a.  Lateral  inflow  or  outflow  is  negligible. 

b.  Channel  width  is  essentially  constant,  i.e.,  9B/3x  » 0 . 

c_.  Energy  losses  from  channel  friction  and  turbulence  are 
described  by  the  Manning  equation. 

d.  Geometry  of  the  section  is  essentially  permanent,  i.e.,  any 
scour  or  fill  is  negligible. 

£.  Bulk  of  the  flood  wave  is  moving  approximately  as  a kine- 
matic wave  which  implies  that  the  energy  slope  is  approxi- 
mately equal  to  the  channel  bottom  slope. 

f_.  Flow  at  the  section  is  controlled  by  the  channel  geometry, 
friction,  and  bottom  slope  and  by  the  shape  of  tne  flood 


8.  If  the  bulk  of  the  flood  wave  moves  approximately  as  a kine- 
matic wave,  the  following  expression  can  be  used  to  transform  the  spa- 
tial gradient  to  time  variations  of  flow  depths: 


where 

y = depth  of  flow,  ft 

x * distance  along  the  river  channel,  ft 
c * kinematic  wave  velocity,  fps 

h = water-surface  elevation  above  mean  sea  level,  ft 
t * time,  sec 

S = effective  bottom  slope  of  the  channel 
o 

r * ratio  of  the  channel  bottom  slope  to  an  average  flood  wave 
slope 

The  last  term  on  the  right  of  the  equation  is  a correction  factor  and 
is  usually  small  for  large  river  applications.  The  kinematic  wave 
velocity  c is  computed  from  the  Kleitz-Seddon  law. 
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where 

B = channel  width  at  the  water  surface,  ft 
Q * discharge,  cfs 

If  the  channel  is  wide  and  prismatic  and  flow  can  be  determined  by 
Manning's  equation,  the  wave  velocity  is 


where  K = 5/3  - ^2A/3B2^iB/dh^ . In  the  above  development,  the  hydrau- 
lic radius  R is  approximated  by  the  hydraulic  depth  D . The  result 
of  the  above  expression  for  wave  velocity  is  that  the  propagation 
speed  of  the  flood  wave  in  the  river  is  a multiple  of  the  river  mean 
velocity.  The  quantity  K varies  from  1.3  to  1.7  for  triangular-  to 
rectangular-shaped  channels,  respectively.  For  a given  channel  geom- 
etry, the  speed  of  wave  propagation  varies  uniquely  with  discharge 
and  stage.  The  wave  speed  is  assumed  the  same  for  the  rising  and 
falling  portion  of  the  input  hydrograph.  Overbank  storage  or  flows, 
backwater  areas,  tributary  effects,  and  other  dynamic  phenomena  are 
not  accounted  for  in  the  model  and  can  become  important  with  complex 
river  systems. 

9.  With  the  definition  of  the  wave  speed  c given  above  and  the 
1* 

method  shown  by  Fread  to  evaluate  r , it  is  possible  to  obtain  the 
following  equation  relating  stage  and  discharge: 


«• 


where 

6h  = change  in  water-surface  elevation  during  the  time  interval 
S At  , in  fps  (6hg  = (h  - h' )/At  , where  h'  is  the  stage  at 
time  t - At). 

Q'  ■ discharge  at  time  t - At  , cfs 
A'  =*  cross-sectional  area  at  time  t - At  , sq  ft 
g * acceleration  due  to  gravity 
At  = small  interval  of  time,  sec 

The  above  equation  is  utilized  in  the  model  as  a basis  for  computing 

discharge  when  the  rate  of  change  of  stage  is  known  or  stage  when  the 

rate  of  change  of  discharge  is  known.  The  unknowns  Q and  h in  the 

above  equation  are  not  expressed  in  an  explicit  manner,  thus  requiring 

a trial -and-correct  solution.  The  solution  of  the  equation  is  by  Newton 

U 

iteration  and  is  explained  fully  by  Fread. 

Model  Implementation 


10.  The  model  (called  DYNMOD  by  Fread)  had  been  implemented  by 
LMVD  on  the  WES  Honeywell  635  computer  by  the  time-sharing  system.  Sev- 
ered. modifications  had  been  introduced  for  ease  of  application.  The 
model  requires  the  input  of  either  a discharge  or  stage  hydrograph  and 
will  compute  the  other  hydrograph.  In  this  application,  stage  hydro- 
graphs were  input  for  calibration  against  field-measured  stage-discharge 
relations  at  the  gaging  stations.  The  model  was  then  used  in  the  dis- 
charge input  mode  to  study  loop  rating  curves  using  several  design 
discharge  hydrographs. 

11.  Channel  input  geometric  data  consist  of:  (a)  the  channel 
effective  bottom  slope  (Sq),  (b)  the  cross-sectional  area  (A)  and  the 
surface  width  (B)  as  a function  of  the  elevation,  and  (c)  the  Manning's 
coefficient  (n)  as  a function  of  the  elevation.  These  are  input  into 
the  model  in  the  form  of  data  tables  at  about  10  elevation  intervals. 

The  area  and  width  were  obtained  from  the  tabulated  1973  and  1975  dis- 
charge observation  data  and  measured  cross  sections  at  the  gaging  sta- 
tions. Plate  1 shows  the  cross  section  at  the  Vicksburg  gaging  site 
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during  the  1973  and  1971*  high-water  hydrographic  surveys.  A low-water 
survey  for  1971  is  also  shown  for  comparison.  Plotted  channel  cross- 
sectional  area  and  widths  for  the  1973  and  1975  flood  years  are  shown  in 
Plate  2 for  the  Vicksburg  gaging  section.  The  graphical  method  of  ex- 
trapolating to  design  stage  levels  is  also  shown.  A preliminary  esti- 
mate of  the  effective  slope  was  obtained  from  the  average  low-water 
plane,  1973  peak  flood  data,  and  channel  bottom  slope  at  the  gaging 
stations.  Trial  computations  required  adjustment  of  this  value  to  give 
a better  representation  of  the  loop  rating  curves.  The  Manning's  n 
values  were  calculated  by  fitting  single-line  rating  curves  at  the 
gaging  stations  using  Manning's  equation  and  known  channel  slope  and 
geometry.  Several  combinations  of  slope  and  resulting  n values  were 
tried. 


PART  III:  MODEL  CALIBRATION  AND  APPLICATION 


Vicksburg.  Miss 


12.  After  model  familiarization,  the  first  step  was  to  recompute 

1* 

a published  test  case  given  in  the  basic  reference  by  Fread.  This  pro- 
vided a check  to  ensure  proper  function  of  the  computer  code.  Several 
improvements  to  the  model  were  developed  during  the  study.  The  data 
input  and  output  structure  was  revised  to  permit  the  use  of  files  for 
ease  of  data  manipulation.  A plotting  capability  using  an  available 
Hewlett-Packard  plotter  was  developed  for  easier  interpretation  of  the 
output  information. 

13.  Most  of  the  model  testing  and  modifications  described  below 
were  made  during  the  calibration  phase  at  the  Vicksburg,  Miss.,  gaging 
station.  Several  combinations  of  river  slopes  and  channel  roughness 
were  tried.  The  river  slope  at  Vicksburg  is  about  0.00006U  ft/ft. 
Computations  were  made  at  effective  slopes  from  O.OOOO6I4  to  0.00001  and 
comparisons  were  made  with  the  measured  stage-discharge  data.  The 
smaller  slopes  tend  to  increase  the  dynamic  loop  width,  but  also  tend 
to  amplify  daily  stage  variations.  A compromise  is  thus  "equired  to 
give  reasonable  model  response  for  reproduction  of  the  1973  rating  curve 
A value  of  0.000032  was  finally  selected.  Table  1 gives  a summary  of 
the  geometric  input  and  channel  roughness  data  at  Vicksburg.  The  data 
tabulated  are  representative  of  the  main  channel  only.  The  plot  in 
Plate  3 shows  the  resulting  computed  and  measured  stage-discharge  curves 
at  Vicksburg;  the  measured  channel  discharge  plotted  has  been  corrected 
for  the  overbank  flow.  It  was  apparent  from  these  results  that  the 
basic  DYNM0D  model  would  not  reproduce  the  magnitude  of  the  loops  as 
compared  with  the  measured  1973  stage-discharge  curve.  The  dynamic  loop 
effect  at  lower  discharges  for  within  bank-full  rising  and  falling 
stages  is  reproduced  very  well  but  the  computed  main  loop  is  only  about 
one  half  the  observed  main  loop. 

lU.  The  model  DYNM0D  was  modified  by  using  different  Manning's  n 
roughness  values  for  rising  and  falling  periods  of  the  flood.  The 


basic  model  uses  a single  set  of  n values  which  along  with  other  input 
data  define  a mean  rating  curve.  The  difference  between  the  computed 
dynamic  stage-discharge  relation  and  the  steady-flow  curve  represents 
the  estimate  of  the  dynamic  terms.  Since  the  results  shown  in  Plate  3 
suggest  that  the  dynamic  terms  are  not  large  enough  to  account  for  the 
observed  loop,  a variable  set  of  n values  was  used  to  obtain  better 
agreement  with  the  observed  1973  flood  data.  In  effect,  the  different 
n values  for  rise  and  fall  were  used  to  account  for  several  factors 
which  cannot  be  properly  evaluated  and  modeled  by  DYNMOD,  such  as 
channel  geometric  changes  from  scour,  changes  in  effective  resistance 
coefficients  due  to  changing  bed  forms , interactions  with  overbank  stor- 
age and  flow,  backwater  storage,  and  inaccuracy  of  estimates  of  the 
dynamic  terms. 

15.  Manning's  n values  to  be  used  for  rising  and  falling  por- 
tions of  the  flood  were  developed  by  drawing  two  rating  curves  through 
the  1973  flood  data  (Plate  3)  and  computing  the  corresponding  n values 
(Table  l).  The  rating  curves  were  extrapolated  to  stages  slightly  above 
the  design  stage  to  provide  roughness  values  for  that  stage  range 
(Plate  5).  After  modifying  the  DYNMOD  model  to  permit  use  of  the  differ- 
ent n values  for  rise  and  fall,  a stage-discharge  curve  was  again 
computed  for  the  1973  flood.  There  was  substantial  improvement  in  agree- 


ment for  the  main  loop  but  the  within  bank- full  loop  was  distorted.  The 


computation  tended  to  switch  between  the  two  sets  of  n values  at  each 
rise  or  fall  of  the  stage  hydrograph.  The  model  was  then  further  modi- 
fied to  allow  the  introduction  of  a switch  to  be  used  to  control  the 
model  to  use  the  rising  or  falling  set  of  n values,  depending  on  pre- 
determined values  of  stage  and  discharge.  This  modification  allows  the 
use  of  the  rising  n values  up  to  a preset  elevation  or  discharge, 
switches  to  the  falling  n values  on  the  first  fall  which  occurs  after 
the  hydrograph  passes  through  the  reference  elevation  or  discharge,  and 
retains  the  use  of  the  falling  n values  for  the  remainder  of  the  com- 
putation even  if  another  rise  occurs.  Using  this  procedure  and  a refer- 
ence elevation  of  9"/  ft  ms  1 (about  50.8  stage  at  Vicksburg),  the  computed 
stage-discharge  curve  for  the  1973  flood  compared  favorably  with  the 
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observed  curve  (Plate  k)  and  the  model  was  considered  calibrated  for 
further  use  at  the  Vicksburg  gaging  station. 

16.  The  design  discharge  at  Vicksburg  for  the  58A-EN  project 
design  flood  (PDF)  is  2,710,000  cfs.  The  peak  flow  during  the  1973 
flood  was  1,960,000  cfs  which  is  about  75  percent  of  the  PDF.  In  order 
to  study  loop  effects  at  design  flows,  an  extrapolation  of  channel  geom- 
etry and  roughness  was  required.  Plate  2 shows  the  method  of  extrapola- 
tion. Channel  area  and  width  were  extrapolated  from  the  1973  and  1975 
data  plots,  taking  into  consideration  a plot  of  the  river  cross  section. 
Overbank  areas  and  widths  were  eliminated  by  vertical  extension  of  the 
main  channel  above  bank-full.  The  areas  and  widths  used  in  the  computa- 
tions above  bank-full  flow  were  representative  of  the  main  channel  only, 
which  was  necessary  because  of  model  limitations. 

17.  An  initial  trial  computation  was  made  by  extrapolating  the 
1973  rating  curve  at  Vicksburg  to  design  flows.  Two  rating  curves 
representative  of  the  rising  and  falling  part  of  the  hydrograph  were 
extrapolated  graphically.  After  additional  study,  it  was  decided  that 
calibrated  channel  Manning's  values  from  the  1973  flood  would  be  extra- 
polated to  design  flows.  Plate  5 shows  a plot  of  channel  roughness 
variation  with  stage  and  the  extrapolation  for  the  rising  and  falling 
part  of  the  hydrograph.  Overbank  flows  during  the  1973  flood  were  sub- 
tracted from  the  total  discharge.  The  PDF  hydrograph  was  also  modified 


to  account  for  estimated  overbank  flows  at  design  discharges.  Overbank 
flow  was  subtracted  from  the  total  design  flow,  based  on  steady-flow 
conveyance  computations  furnished  by  the  Vicksburg  District.  The  method 
of  flow  partition  was  based  on  computations  made  using  the  steady-flow 
backwater  computations  using  HEC-2.  The  peak  flow  at  Vicksburg,  for 
example,  was  reduced  from  2.71  million  cfs  to  2.56  million  cfs.  The 
resulting  design  hydrographs  are  tabulated  in  Table  2. 

18.  After  the  calibration  and  modifications  previously  discussed 
had  been  made  to  DYNMOD,  two  computations  were  made  using  as  input  the 
project  design  flow  hydrograph  58A-EN  and  a composite  hydrograph  which 
follows  the  1973  flood  through  its  peak  flow  of  about  1.96  million  cfs, 
then  assumes  a subsequent  rise  to  the  design  flow  of  2.56  million  cfs 


following  the  shape  of  the  58A-EN  hypothetical  flood  (Table  2).  Re- 
sults of  the  computations  are  shown  in  Plate  6.  The  double-valued 
Manning's  n computation  was  used  to  simulate  the  effects  of  a 1973 
type  of  flood  followed  by  the  design  flood,  using  the  falling  n values 
for  flows  above  the  1973  peak  and  up  to  design  discharges.  The  run 
using  only  the  rising  n values  and  the  58A-EN  flood  simulates  the 
Vicksburg  rating  curve  if  the  channel  roughness  had  not  increased  for 
corresponding  stages  during  the  design  flood.  The  difference  between 
the  two  curves  at  design  flows  is  a measure  of  the  dynamic  loop  effect 
(plus  other  unaccounted  for  factors  such  as  overbank  effects  and  bed- 
form  changes).  Assuming  that  the  higher  channel  roughness  would  be 
appropriate  if  the  design  hydrograph  followed  the  1973  peak  flows,  the 
difference  in  the  two  curves  can  be  used  for  estimating  loop  effect. 

The  computations  indicate  a difference  of  about  1.7  ft  in  the  peak 
stages  with  and  without  the  loop  effect  at  Vicksburg. 

Helena,  Ark. 

19.  A plot  of  the  Helena  discharge  range  cross  section  is  shown 
in  Plate  7 for  the  1973  and  1975  high-water  surveys.  The  1973  cross 
section  was  used  to  calculate  the  cross-sectional  area  and  channel  width 
and  these  values  are  plotted  in  Plate  8.  Other  width  and  area  measure- 
ments made  in  1973  and  1975  flood  years  at  the  discharge  range  are  also 
plotted  in  Plate  8. 

20.  A calibration  procedure  was  used  at  Helena  similar  to  that 
developed  for  the  Vicksburg  gaging  station.  The  final  results  of  sev- 
eral trial  computations  at  Helena  are  shown  in  Plate  9.  The  switch  from 
rising  to  falling  n values  was  set  to  U8.6-ft  stage  and  1,550,000  cfs 
for  Helena  to  provide  the  best  fit  to  the  1973  flood  discharge  data. 

The  plot  indicates  good  agreement  with  the  channel  discharge  measure- 
ments. The  channel  geometric  data  and  Manning's  n computations  are 
presented  in  Table  3.  The  cross-sectional  area  and  width  were  based  on 
the  average  data  plotted  in  Plate  8.  Channel  area  and  width  at  design 
stage  were  obtained  by  plot  extrapolation.  The  calculated  set  of 


Manning's  n is  shown  plotted  in  Plate  10.  Several  channel  slopes  were 
tried  from  the  apparent  channel  slope  of  0.00011  down  to  a value  of 
0.000032  ft/ft;  a slope  of  0.000064  was  finally  selected.  The  effective 
channel  slope  and  computed  roughness  values  at  Helena  are  about  twice 
those  at  the  Vicksburg  gaging  station.  The  channel  roughness  values  at 
Helena  appear  to  be  in  agreement  with  other  recent  studies'* of  Missis- 
sippi River  channel  roughness.  The  Vicksburg  values  seem  to  be  low. 
Increasing  the  slope  at  Vicksburg  to  0.000064  would  increase  n values 
by  *^2  or  4l  percent,  which  would  still  be  lower  than  the  Helena  data. 
Several  factors  not  accounted  for  in  the  model,  such  as  Yazoo  backwater 
area  effects  and  changing  river  sinuosity  with  stage,  could  be  the  rea- 
son for  these  apparent  discrepancies. 

21.  Computations  were  run  using  the  58A-EN  design  flood  at  Helena 
and  a design  discharge  hydrograph  made  up  of  1973  daily  Helena  dis- 
charges phased  into  the  58A-EN  hypothetical  flood  at  Helena.  These  data 
are  tabulated  in  Table  4.  The  design  hydrograph  was  corrected  for  over- 
bank flows  based  on  backwater  computations  by  the  Memphis  District.  One 
run  was  made  using  the  1973/58A-EN  hydrograph  with  different  Manning's 

n values  on  the  rising  and  falling  parts  of  the  hydrographs.  Plate  11 
presents  the  results.  Similar  reasoning  at  Helena  to  that  developed 
at  Vicksburg  relative  to  the  sequence  of  flood  flows  and  channel  rough- 
ness suggests  that  the  difference  between  the  two  curves  at  design  flows 
is  a measure  of  the  dynamic  loop  and  other  effects.  The  difference  in 
the  two  curves  can  be  used  for  estimating  loop  effect  and  indicates  the 
total  loop  effect  to  be  about  2.0  ft  at  Helena. 

Baton  Rouge,  La. 

22.  The  channel  cross  section  at  the  Baton  Rouge  gaging  range  is 
plotted  in  Plate  12  based  on  the  January  1974  hydrographic  survey. 

There  is  no  overbank  area  at  Baton  Rouge.  The  calibration  procedure  was 
modified  at  the  Baton  Rouge  gaging  station.  A study  of  the  1973  data 
indicated  large  scatter  in  the  channel  cross-sectional  area  and  width 
with  stage  (see  Plate  13).  During  the  1973  flood,  discharge  measurements 
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were  made  only  during  rising  stages  from  39  ft  to  42  ft  and  falling 
stages  to  about  32  ft.  The  loop  rating  curve  is  very  hard  to  define  and 
extrapolate  to  design  flows  without  discharge  data  to  show  the  rising 
and  falling  hydrographs  at  a wide  range  of  river  stages.  Additional 
data  were  reviewed  from  the  1974  and  1975  high  water  but  these  showed 
similar  data  scatter  and  a narrow  range  of  river  stage.  In  addition, 
measurement  techniques  and  discharge  range  location  apparently  were 
varied  during  1974  and  1975. 

23.  Area  and  width  data  from  the  1950  flood  were  also  plotted  in 
Plate  13  and  showed  much  less  scatter.  Discharge  measurements  were 
made  over  a wider  range  of  stages  during  the  1950  flood.  The  1950  dis- 
charge range  was  located  3.1  miles  downstream  from  the  1973  discharge 
range  which  was  at  the  U.  S.  Highway  190  Bridge  or  river  mile  230.8, 

1962  AHP.  It  was  decided  to  calibrate  the  model  using  both  1950  and 
1973  flood  years.  The  channel  data  for  the  1950  discharge  measurements 
appeared  adequate.  A study  of  the  1973  discharge  data  indicated  that 
most  of  the  measurements  were  made  at  a discharge  range  about  200  ft 
downstream  of  U.  S.  Highway  190  Bridge.  The  channel  cross  section  shown 
in  Plate  12  was  used  as  the  best  source  of  geometric  information  for  the 
1973  flood.  Table  5 presents  the  adopted  geometric  data  for  the  1950 
and  1973  floods.  Note  that  Manning's  n values  were  not  varied  for 
the  rising  and  falling  portion  of  the  hydrograph.  As  suggested  above, 
the  basic  data  did  not  warrant  the  use  of  multivalued  Manning's  n co- 
efficients. As  had  been  done  previously  at  Vicksburg  and  Helena  gaging 
stations,  several  channel  slopes  were  tried  from  0.00001  to  0.000032 
ft/ft.  The  value  of  0.000016  ft/ft  was  used. 

24.  The  calibrations  of  the  model  to  the  1950  and  1973  measured 
discharge  data  are  shown  in  Plates  14  and  15,  respectively.  The  1950 
calibration  gave  poor  reproduction  of  the  rising  part  of  the  hydrograph 
up  to  about  the  4l-ft  stage.  Agreement  was  adequate  at  higher  stages. 
The  1973  calibration  showed  large  fluctuations  and  very  complex  stage 
discharge  relations  at  stages  above  39  ft.  It  is  speculated  that  gate 
operation  at  the  Morganza  Floodway  caused  these  variations.  The  result- 
ing Manning's  n values  are  plotted  in  Plate  16. 


25.  The  project  design  flood  hydrograph  58A-EN  at  Baton  Rouge  was 
furnished  by  the  New  Orleans  District  and  is  shown  in  Table  6.  Because 
of  the  designed  operation  of  the  Morganza  Floodway  during  a large  flood, 
flows  at  Baton  Rouge  would  be  kept  below  1,500,000  cfs.  In  order  to  pro- 
vide a realistic  hydrograph  that  gives  some  weight  to  flow  sequences,  a 
composite  design  hydrograph  was  obtained  at  Baton  Rouge.  This  hydro- 
graph was  calculated  by  combining  the  58A-EN  hydrograph  with  increased 
1973  daily  flows  from  April  18  through  June  30.  The  daily  flows  at 
Baton  Rouge  were  multiplied  by  the  ratio  of  1,500,000  to  1,31*7,000 
(l.llU)  assuming  the  design  flow  maximum  of  1,500,000  cfs  occurred  on 

18  April  1973.  This  was  the  first  discharge  measurement  available  at 
Baton  Rouge  after  the  opening  of  the  Morganza  Floodway  on  17  April  1973. 
This  hydrograph  is  an  attempt  to  simulate  the  river  behavior  at  Baton 
Rouge  which  is  representative  of  1973  conditions  at  the  design  flow  of 
1,500,000  cfs.  Both  design  hydrographs  are  shown  in  Table  6. 

26.  Computations  were  made  with  both  design  hydrographs  and  are 
shown  in  Plates  17  and  18.  The  design  flood  58A-EN  computations  show  a 
dynamic  loop  effect  of  about  0.7  ft  near  the  peak  flow  and  a maximum 
stage  of  1*5  ft.  The  composite  1973/58A-EN  design  flood  computations 
result  in  a loop  effect  somewhat  less  than  1.0  ft  and  a maximum  stage 
of  1*5.2  ft.  Results  of  both  computations  indicate  a dynamic  loop  effect 
of  about  1.0  ft  at  Baton  Rouge.  The  amount  of  dynamic  effect  in  the 
Mississippi  River  in  the  Old  River  to  Baton  Rouge  reach  during  a design 
flood  will  be  greatly  dependent  on  the  mode  of  operation  of  the  Old 
River  Control  Structure  and  Morganza  Floodway.  If  these  operations  are 
expected  to  be  greatly  different  than  design  or  1973  conditions,  dynamic 
effects  from  gate  manipulations  may  cause  higher  stages.  A more 
thorough  analysis  using  a complete  dynamic  model  would  be  needed  to 
ascertain  stage  response  to  variable  gate  structure  operations. 
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PART  IV:  DISCUSSION  OF  RESULTS 
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27.  For  high  flows  at  the  Vicksburg  and  Helena  gaging  stations, 
the  computation  follows  the  extrapolated  n values  for  either  the  ris- 
ing or  falling  portions  of  the  1973  flood  (Plates  6 and  ll).  The  differ- 
ence between  computations  for  the  different  flood  sequences  is  attributed 
to  loop  effect.  The  rationale  supporting  this  type  of  projection  is: 
had  the  last  rise  to  design  flow  occurred  in  1973  or  a similar  antece- 
dent flood,  the  roughness  values  compatible  with  the  falling  portion  of 
the  antecedent  flood  would  be  more  appropriate  during  that  rise  than  the 
rising  n values  associated  with  the  lower  channel  roughness.  The  use 
of  the  modified  Fread  model  DYNMOD  in  this  fashion  is  to  provide  for  a 
judgment  prediction  of  the  loop  effect  rather  than  as  a precise  predic- 
tive or  modeling  tool.  Nevertheless,  the  loop  adjustment  suggested 

by  these  projections  is  believed  to  be  valuable  in  supplementing  the 
steady-state  flow  line  computed  from  HEC-2.  In  other  words,  the  dynamic 
model  was  used  in  the  analogous  manner  as  steady-state  computations  were 
made  and  the  model  is  believed  consistent  and  comparable  in  accuracy. 

28.  The  setting  of  the  project  design  flow  line  is  to  be  used  for 
the  design  of  levee  crest  heights  and  other  water-control  structures 
along  the  main  stem  of  the  Mississippi  River.  Previous  flow  line  compu- 
tations in  1957  had  been  made  using  steady-state  backwater  computations 
and  maximum  design  discharges.  The  rating  curves  used  were  a series  of 
single-valued  curves  based  on  field  measurements  during  the  19^5  and 
1950  floods.  This  type  of  computation  was  repeated  using  the  1973 
flood  data  as  a part  of  the  flow  line  reanalysis.  The  previous  LMVD 
report  recognized  the  need  to  allow  for  dynamic  effects  over  and  above 
the  steady-state  flow  line.  Steady-state  (or  equivalently,  single  line) 
rating  curves  are  usually  an  average  between  the  rising  and  falling  part 
of  the  hydrographs.  A loop  magnitude  allowance  that  would  be  consistent 
with  the  steady-state  flow  line  would  then  be  about  one  half  the  total 
dynamic  loop  magnitude. 

29.  Results  obtained  with  the  DYNMOD  model  suggest  that  channel 
dynamic  effects  are  small  near  design  flow  conditions  even  though  they 
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may  be  substantial  at  smaller  flows.  Even  at  the  lower  flows,  dynamic 
effects  do  not  appear  to  account  for  the  entire  loop  observed  during  the 
1973  flood.  With  the  model  used  herein,  it  was  necessary  to  lump  the 
effects  of  several  factors  into  the  roughness  coefficient.  This  empha- 
sizes the  need  for  accuracy  in  extrapolating  roughness  coefficients  to 
design  flows  and  for  investigations  to  evaluate  the  effects  of  other 
complex  river  factors.  Better  definition  of  the  effects  of  some  of 

these  factors  can  be  obtained  with  more  detailed  studies  using  modeling 

7 

methodologies  such  as  developed  for  the  junction  near  Cairo,  111. 

Fread  has  recently  suggested  a methodology  for  more  accurate  treatment 

g 

of  the  overbank  effect  that  includes  overbank  storage  and  flow  between 
main  channel  meanders.  Evaluating  the  effects  of  changing  channel 
geometry  and  roughness  due  to  bed  changes  will  be  more  difficult. 

Because  of  the  difficulty  in  predicting  bed  changes  and  their  effects, 
it  would  appear  advisable  to  use  available  modeling  methodology  to 
accurately  evaluate  all  factors  associated  with  a fixed  geometry  system 
before  concluding  that  the  sediment  transport  phenomena  must  be  treated 
to  obtain  predictions  of  acceptable  accuracy. 
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PART  V:  CONCLUSIONS 


30.  As  a result  of  the  investigation  reported,  the  following  con- 
clusions have  been  formulated: 

a.  The  study  of  loop  rating  curves  on  the  Mississippi  River 
has  shown  that  the  Fread  model  DYNMOD  will  account  for  ' 
about  one  half  of  the  total  loop  magnitude  when  using  only 
a single  set  of  n values. 

b.  With  computations  using  a single  set  of  n values,  peak 
design  discharges  produced  the  same  maximum  stage  regard- 
less of  the  shape  of  the  discharge  hydrograph,  indicating 
that  the  dynamic  factors  alone  do  not  properly  account 
for  total  loop  effect. 

c_.  The  sequence  of  hydrograph  flows  does  not  have  an  appre- 
ciable effect  on  the  rating  curve  predicted  by  DYNMOD 
using  a single  set  of  n values.  The  dynamic  loop 
effect  on  the  rating  curve  is  small  at  and  near  design 
flows  where  rates  of  rise  are  usually  small. 

d.  The  introduction  of  double-valued  Manning's  n for  rising 
and  falling  hydrographs  with  a preset  switch  is  necessary 
for  simulating  the  1973  flood  event  or  other  multiple 
peaked  flood  events. 

e_.  Using  the  modified  model  DYNMOD,  it  is  possible  to  simu- 
late dynamic  effects  at  design  flows,  assuming  that  the 
loop  behavior  during  the  1973  flood  could  have  occurred 
at  design  flows. 

f_.  The  rating  curves  computed  assuming  single- valued  rough- 
ness with  stage  when  compared  with  the  rating  curves  com- 
puted assuming  double-valued  roughness  give  a reasonable 
extrapolation  of  the  full  dynamic  loop  effect.  Several 
hydraulic  phenomena  such  as  tributary  flows , channel- 
overbank  interaction,  and  channel  bed  changes  are  lumped 
into  the  roughness  values  used  in  the  model  to  simulate 
the  1973  flood. 

£.  Considering  the  method  in  which  the  steady-state  model 

calibration  was  accomplished,  and  the  computational  proce- 
dure used  in  the  modified  model  DYNMOD,  it  is  concluded 
that  one  half  of  the  total  loop  effect  should  be  used  as 
the  allowance  for  loop  effect  in  establishing  levee  grades 
throughout  the  studied  reach.  Averaging  the  values  cal- 
culated at  Vicksburg  and  Helena,  and  using  one  half  of 
the  average,  the  total  loop  effect  to  be  added  to  steady- 
state  stages  would  be  1.0  ft. 
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h.  An  allowance  of  approximately  1.0  ft  should  be  added 
the  flow  line  from  about  Helena  to  Red  River  Landing 
amount  of  0.5  ft  should  be  sufficient  from  Red  River 
ing  to  about  Baton  Rouge. 
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Input  Channel  Data,  Vicksburg,  Miss. 
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Table  3 

Input  Channel  Data,  Helena,  Ark. 


Bottom  Slope  * 0.00006U  ft /ft 
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Table  5 

Input  Channel  Data,  Baton  Rouge.  La. 
Bottom  Slope  = 0.000016  ft /ft 


Stage 

ft 

Elevation 
ft  msl 

Area 
x 1000 
sq  ft 

Width 

ft 

Manning's 

n 

1950  Flood  Data 

20 

19.1*8 

1U0 

3090 

0.0201 

25 

2U. 1*8 

155 

3110 

0.0192 

30 

29.1*8 

170 

3130 

0.0183 

32.5 

31.98 

178 

3180 

0.0180 

35 

31*.  1*8 

185 

3360 

0.0168 

1*0 

39.1*8 

201 

3650 

0.011*3 

1*5 

1*1*.  1*8 

220 

3690 

0.0133 

1*7 

1*6. 1*8 

226 

3690 

0.0129 

1973  Flood  Data 

20 

19.1*8 

139 

2330 

0.021*0 

25 

2l*.  1*8 

152 

2380 

0.0222 

30 

29.1*8 

161* 

21*90 

0.0201 

35 

3U.U8 

177 

2820 

0.0175 

1*0 

39.1*8 

190 

2880 

0.0152 

1*5 

1*1*. 1*8 

205 

2920 

0.0138 

1*7 

1*6. 1*8 

211 

2930 

0.0131* 

I 


I 


Table  6 

Design  Flood  Hydrographs.  Baton  Rouge,  La. 
Daily  Discharge,  1000  cfs 
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